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Weak lensing peak count as a probe of f{R) theories 
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ABSTRACT 

Weak gravitational lensing by galaxy clusters on faint higher redshift galaxies 
has been traditionally used to study the cluster mass distribution and as a tool to 
identify clusters as peaks in the shear maps. However, it becomes soon clear that 
peaks statistics can also be used as a way to constrain the underlying cosmological 
model due to its dependence on both the cosmic expansion rate and the growth rate of 
structures. This feature makes peak statistics particularly interesting from the point 
of view of discriminating between General Relativity and modified gravity. Here we 
consider a general class of f{R) theories and compute the observable mass function 
based on the aperture mass statistics. We complement our theoretical analysis with 
a Fisher matrix forecast of the constraints that an Euclid - like survey can impose 
on the /(i?) model parameters. We show that peak statistics alone can in principle 
discriminate between General Relativity and f{R) models and strongly constrain the 
/(i?) parameters that are sensitive to the non-linear growth of structure. However, 
further analysis is needed in order to include possible selection function in the peaks 
redshift determination. 
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1 INTRODUCTION 

The gravitational field of galaxy clusters distorts the im- 
ages of background faint galaxies possibly leading to the 
formation of specta cular giant arcs (as first obser v ed by 
ISoucail et all Il987l '). As pointed out by IWebsterl (|l985l ) 
in a pioneering work, the most likely effect is, however, a 
variation in the ellipticity distribution of the background 
galaxies which can then be used to recons truct the mass dis- 
tribu ti on of the lens i ng clu ster (see, e.g.. iKaiser fc Squires I 
1 19931 : H aiser et al. I I1995I '). Such a mass reconstruction 
method has now become quite popular JClowe et al. 19981: 
Radovich et al. I I2OO8I : iRomano et al. I I2OI0I : iHuang et al. I 



be used as a way to find dark mass concentrations thus 
offering the opportunity of assembling a mass selected 
catalog of haloes. In particular, the aperture mass statistics 
(|Schneider J 1996|) has emerged as a v aluable way to find 
clust e rs (IHetterscheidt et al. I l2005l : iGavazzi fc Soucail I 



I2OO7I : ISchirmer et al. 



20 111 ) and is actually the most efficient one for intermedi- 



ate and high redshift clusters where dynamical methods fail 
because of the difficulties in measuring galaxy redshits. 

The search for coherent image alignments can also 



I2OO7I ) and measuring their mass 
function (|Dahle 1 12006| ) . On the other hand, although not 
originally conceived as a cosmological tool, the aperture 
mass statistics allows to severely constrain cosmological 
models. Indeed, the number counts of peaks in the weak 
lensing maps is determined by both the cosmic expansion 
rate and the growth rate of structures (entering through 
its effect on the theoretical mass function). Peaks statistics 
has therefore emerged as a promising t ool to discriminate 
among different dark energ y models (iBartelmann et al. I 



20021: iMarian et al 



2001 



Maturi et al. II2OIOI : IKratochvil et al 



Dietrich fc Hartlap I I2OI0I: 



primordial non - Gaussianity ( Maturi et al. |[201ll ) 



2OIOI) and constrains 
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Actually, dark energy is not the only way to fit 
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the wide amount of data indicating the present day uni- 
verse as spatiaUy flat and undergoing accelerated expan- 
sion (jWeinberg et al. II2OI2I ). Indeed, rather than represent- 
ing the evidence of a missing source in the cosmic energy 
budget, cosmic speed up may also be read as the first signal 
of breakdown of our understanding of the laws of gravity. 
Modified gravity theories has therefore recently attracted a 
lot of interest and have been shown to be able to fit the data 
with the same accuracy as most dark energy models. In this 
framework, f{R) theories represent one of the most natu- 
ral generalization of General Relativity, the basic idea being 
to replace the scalar curvature R in the gravity Lagrangian 
with a generic function f{R)- Dating back to Eddington, this 
idea has been first reconsidered in the '80s thanks to the pos- 
sibility of recovering inflation without any scalar field and 
has then found a renewed interest after the discovery of the 
acceleration for its ability to achieve cosmic speed up in a 
matter only universe (see. e.g., | de Felice fc Tsuiikawa ' 



j i m a 
2OIOI : 



ICapozziello fc de Laurentis 1 l201lf and refs. therein). It is ac- 
tually possible to show that, for any dark energy model, it 
is possible to work out an f{R) counterpart providing the 
sam e background expansion, i.e. the same H(z), as the g iven 
one (|CapozzieUo et al. ll2005l : lMultamaki fc Vilia Il2006l 'l. On 
the contrary, these equivalent models can be discriminated 
because of the different growth of perturbations driven by a 
scale and redshift dependent effective gravitational constant 
and the non vanishing difference between the two Bardeen 
potentials. 

As already reminded above, the number counts of weak 
lensing peaks is indeed sensible to both the expansion rate 
and the growth of structures so that it comes out as an ideal 
tool to discriminate among dark energy and f{R) models 
sharing the same H{z) expression. Here, we will therefore 
investigate whether this is indeed the case for a popular 
f{R) model which has already been shown to excellently 
fit a wide astrophysical dataset. Such a preliminary study 
can highlight the power of peaks statistics in discriminating 
among dark energy and modified gravity, but it is not imme- 
diately related to what one could actually infer from obser- 
vations. We will therefore rely on the Fisher matrix analysis 
to forecast the accuracy a realistic survey can achieve on the 
f{R) model parameters thus filling (although in a simplified 
way) the gap between theory and observations. 

The plan of the paper is as follows. In Sect. 2 we review 
the basic of f{R) theories and how one can compute the 
mass function taking care of the peculiarities of this modified 
gravity model. Sect. 3 is then devoted to the mass aperture 
statistics discussing all the steps needed to compute both the 
signal and the noise in weak lensing maps. The expected 
number of counts of peaks in f{R) theories and how this 
depends on the model parameters are presented in Sect. 4, 
while the Fisher matrix forecasts are given in Sect. 5. Sect. 6 
summarizes the main results and our conclusions. 



in the literature. We here first briefly review the basics of 
f{R) theories and then explain the method used to evaluate 
the theoretical mass function. 



2.1 f{R) cosmology 

In the framework of the metric approach, the field equations 
are obtained by varying the gravity action 



S = d xy^—g 



fiR) 



(1) 



(2) 



16-kG 

with respect to the metric components. We obtain 

where R is the scalar curvature. Cm is the standard matter 
Lagrangian with Tf^v the matter stress - energy tensor, and 
the prime (the dot) denotes derivative with respect to R 
(time t). For f{R) = R — 2A, one obtains the usual Einstein 
equations with a cosmological constant A, while, in the gen- 
eral case, a further scalar degree of freedom is introduced. 

In a spatially fiat homogenous and isotropic universe, 
some convenient algebra allows to rearrange the field equa- 
tions in such a way that a single equation for the Hubble 
parameter H = a/ a is obtained. Assuming dust as grav- 
ity source and introducing E = H{z)/Ho, it is then only a 
matter of algebra to get : 



^2^^^ ^ [{l + zf + {U'-m^f)/6] 



f -m'2{l + z){d^/dz) 



2 J.// rf^C , 4 r/rr { d(, 



dlnE(z) 



E^z) 



(1 + -)- 



dln(l + z) 



2mV - W 



3(H-z)2 



1 + z dz 
(4) 



where z = 1/a — 1 is the redshift, ^ = R/rn^ and 



2 _ (87rG)VM(z = 0) 



(8315 Mpc)" 



0.13 



(5) 



is a convenient curvature scale which depends on the present 
day values of the matter density parameter Q,m and the 
Hubble constant h = i/o/(100 km/s/Mpc). 

By inserting Eq.Q into Eq.(|4]), we get a single sec- 
ond order nonlinear differential equation for ^(z) that can 
be solved numerically provided f{K) and the initial condi- 
tions are given. These latter can be conveniently expressed 
in terms of the present day values of the deceleration [q = 
—H^^d/a) and jerk {j = H~^ d /a) parameters. To this end, 
we first remember that, in a FRW spatially flat universe, the 
curvature scalar R reads : 



2 FOURTH ORDER GRAVITY 

Being a straightforward generalization of Einstein GR, 
fourth order gravity theories have been investigated almost 
as soon as the original Einstein theory appeared. It is, there- 
fore, not surprising that the corresponding field equations 
and the resulting cosmology have been so widely discussed 



R = 6{H + 2H'^) 

so th at we get for the present-day values (|Capozziello et alTI 
l2008t ) : 



Ro = 6H!{l-qo) , Ro = 6H^{jo-qo-2) 
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It is then only a matter of algebra to get the initial conditions 
for 5 as 



C(z = 0) = (6/!^M)(l-go) 
d^/dz{z = 0) = (6/nM)(jo - go - 2) 



(6) 



Because of the definition of ^, Eq.Q is a single fourth - order 
nonlinear differential equation for the scale factor a{t) so 
that we need to know the values of the derivatives up to the 
third order to determine the evolution of a{t) thus explaining 
why the jerk parameter also enters as a model parameter. 

A key role in fourth order theories is obviously played 
by the functional expression adopted for f{R)- We choose 
here the Hu & Sawicki (2007, hereafter HS) model setting : 



f{R) = R-m 



2 ci(E/m^)" 

l + C2(i?/m2)" 



(7) 



where is given by ((Sj, and (n, ci, C2) are positive dimen- 
sionless constants. Note that, since 



lim f{R) ~ i? - —m" + 



C2 



R 



we recover an effective A term in high curvature {m^ /R — > 0) 
environments. In particular, in the limit R >> m^, the 
expansion rate H, in the early universe, will be the same 
as in ACDM with an effective matter density parameter 
^M,eff = 6c2/(ci + 6C2) that guarantees that the nucle- 
osynthesis constraints are satisfied. 

The HS model is determined by the three parameters 
{n,ci,C2), but it is actually convenient to replace (ci,C2) 
with (qo,.?o). To this end, we follow the method detailed in 
ICardone et aTl (|2012l ') and introduce the further parameter 
£ = log (/ijo — 1) = log [f'{z = 0) — 1] which quantifies the 
present day deviation of f{R) from GR. Note that, since we 
expect /_Ro — 1 << 1, e will typically be negative and large. 

Inserting the HS f{R) into Eqs.© and (Q and setting 
the model parameters, one can solve the fourth order non- 
linear differential equation for the scale factor a{t) and then 
work out the dimensionless Hubble rate E{z). In order to 
speed up the computation, however, we use here the follow- 
ing approximated solution jCardone et al. ||2012| ') : 

r £{z)Ecpl(z) + [1-£{z)]Ea{z) ziizA 
E{z) = <^ (8) 



where 



E'r 



= nAf(l+2)' + (l-nA/)(l + 2)'^'+™''+""'^e-^ (9) 



is the Hubble rate for the phenom enological CPL 
l|Chevallier fc Polarski lliooH : iLinder |[20o3 ) parametrization 
of the effective dark energy fiuid equation of state, namely 
■w{z) — wq + Wail — a), while Ea{z) — Ecpl{z,wo = 
— l,Wa — 0). Finally 



£{z) ^^ei{z ■ 



ZAj 



(10) 



is an interpolating function with ei and za depending on 
(Sli\/, go, jo, 71, e). This approximating function excellently 
reproduces the numerical solution whatever the model pa- 
rameters {QM,<lo,jo,n,e) are with a rms error which is far 
lower than 0.1% over the redshift range (0, 1000). A subtle 
remark is in order here concerning the value of Qm- Indeed, 
while for z ^ za, we use the actual matter density parame- 
ter, the effective one must be used for z ^ za- Therefore, a 
discontinuity in za is formally present in our approximation. 
Actually, it is easy to show that, for all reasonable model pa- 
rameters, Qm and Q.M,eff are almost perfectly equal so that 
the discontinuity can not be detected at all and E{z) is, to 
all extents, a continuous function. 



2.2 Mass function for f{R) gravity 

The extra scalar degree of freedom introduced by the mod- 
ified gravity Lagrangian mediates an enhanced gravita- 
tional force on scales smaller than its Compton wavelength. 
In order to hide this boost from local tests of gravity, 
viable .f(R) models re s ort to the chameleon mec hanism 
iMota fc Barrow 1 12004 iKhourv fc Wehman I [200I) which 
makes the field Compton wavelength depending on the en- 
vironment gravitational potential. As a consequence, for vi- 
able f{R) models, the gravitational potential reduces to the 
Newtonian one on Solar System scales, while nonlinearities 
in the field equations appear for galaxy scale systems. From 
the point of view of the mass function, the enhanced gravi- 
tational force has a particularly strong impact on the abun- 
dance on intermediate mass haloes. Indeed, on one hand, 
the extra force increases the merging rate of low mass haloes 
into intermediate ones, while the chameleon shuts down the 
merging of these latter into highest mass ones. 

An efficient method to take the chameleon effect into 
account and work out a mass function in agreement with 
simulations has been recently proposed bv iLi fc Hu I (|201ll ) 
based on a simple interpolation scheme. Following their ap- 
proach, we model the mass function as : 



A/'(ln M„i 



Pm{z = 0) dlnz^ 



(11) 



Mvir din M^i 

where ip{i') is the number of halos with the ratio between 
density contrast and variance of the perturbations equal to 
v = Sc/a-{Myir), Sc is the critical overdensity for collapse 
and a the variance of the perturbations on the scale R cor- 
responding to the mass Mvir, given by 



a^[RiM,,r)] = 



1 



(27r)3 



Ps{k)\W{kR)\'^d^k 



(12) 



Here, W{kR) is the Fourier transform of the spherical top 
hat function, while the density power spectrum Ps{k) may 
be computed as : 



Ps{k,z) ^ Ak"^-''T''{k)D\k,z) 



(13) 



with T{k) the transfer function, D{k,z) — 5{k, z) / 5{k,Q) 
the growth factor normalized to 1 at present day and A a 
normalization constant conveniently expressed as function 
of CTs, the present day variance on the scale R — 8h~^ Mpc. 
Since f{R) models reduce to GR in the prerecombination 
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Figure 1. Effective gravitational constant for the fiducial HS model as a function of z (left panel) for k = 0.05 (blue), 0.10 (red), 0.15 
(black) and k (right panel) for 2 = (blue), 0.5 (red), 1.0 (black). 



epoch, we can use for T (k) the analytical approximation in 
lEisenstein fc HtTI (| 19981 ) although this latter has been ob- 
tained assuming GR validity. On the contrary, the growth 
factor for /(J?) models gains a scale dependence which is 
not present for GR dark energy models. Indeed, the evolu- 
tion of the density perturbations, 5 = 5pm/pm, in fourth 
order gravity models is ruled by 



5 + 2H5 - ATig^ff{a, k)pM5 = 
with k the wavenumber and 



Gess{a,k) 



G l + 4{k^/a')[f"{R)/f'(R)] 
f'{R)l + 3{kya^)[f"{R)/f'{R)] 



(14) 



(15) 



the scale depe ndent effective gravitational constant 
l|Tsuiikawa 1120071 ). Note that, for f{R) = R-2K, g^j; re- 
duces to G and the growth factor is no more scale depen- 
dent. As shown in Fig. [3 setting the model parameters to the 
fiducial values, Gef / only slightly deviates from the Newton 
constant. Nevertheless, this deviation is both redsfhift and 
scale dependent so that it has an impact on the final matter 
power spectrum. 

In order to compute the mass function through Eq.([TT|, 
one has to choose an expression for vtpiv). We again fol - 
low iLi fc Hul (|201ll ) and adopt the lSheth fc Tormen I (|l999l ) 
function setting 



[l-h(ai^")-^]exp(-aj/72) 



(16) 



with {A,a,p) = (0.322,0.75,0.3). It is, finally, worth stress- 
ing that both 5c and A„ir (which enters the conversion from 
Mmr to Rvir) actually depends on the matter density pa- 
rameter Q.M and the physics ruling the collapse and virializa- 
tion of perturbations. Previous attempts to fit the numeri- 
cally derived mass function searched for a /(i?) derivation of 
(5c, A„ir) and then use the mass variance a^M^ir) evaluated 
from the f{R) linear pow er spectrum. On the contrary, here 
we follow the approach of lLi fc Hu I ()201ll ) holding (tJc, A^ir) 
fixed to the values predicted for the GR model having the 
same background expansion as the /(i?) one. We then use 



for the mass variance an interpolation between the f{R) and 
GR ones, i.e. a{A'Ivir) is estimated as: 



""^^-^^ l-f(M/M..)^ ^''^ 

with apoG and acR computed from the f{R) and GR linear 
power spectra and {a,Mth) two interpolation parameters 
given by 



r a = 2.448 



Mth = 1-345 X 10' 



13 / l/flo 



10- 



3/2 



h'^ Me 



(18) 



Finally, we follow [Bryan fc Norman] (|l998l ) to get {5c, ^vir) 
as a function of Qm and the dimensionless Hubble rate E{z). 

2.3 Fiducial model and theoretical mass function 



The a pproximated formula ((HJ for E{z) and the iLi fc Hul 
(|201ll ) formalism allows us to estimate the mass function 
for the HS f{R) model provided we set the parameters 
(JIa/, h, go, Jo, n, e). Moreover, one also has to set the slope of 
the power spectrum through nps and fix its normalization 
through (jg. Since we are mainly interested in how the mass 
function looks in f{R) models and whether it is possible to 
discriminate among the HS model and its corresponding GR 
based dark energy model with the same expansion rate, we 
will not vary all the parameters at play. On the contrary, we 
will first fix the expansion rate to that of a GPL model with 

(fijv/, h, wo, wa) = (0.273, 0.703, -0.95, 0.0) . 

Second, we fix the shape and the amplitude of the mat- 
ter power spect rum setting (np S jO's) = (0.966,0.809) in 
agreement with iKomatsu et al. I (|201lf ) constraints. As a 
fiducial HS model, we choose the one obtained by setting 
{n,e) — (1.5, —6.0) since it best mimics the GPL expansion 
rate. Moreov er, this model is quite similar to the best fit 
one found bv lGardone et ai~l (|2012l ) fitting a large dataset 
including SNela, GRBs, BAOs and WMAP7 distance priors. 
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Figure 2. Left. Mass function (MF) for the fiducial HS model for z = (blue), z = 0.5 (purple), z = 1.0 (red) shifting by 0.25 (0.50) 
dcx the MF at z = 0.5 (z = 1.0) for a better visualization. Right, Pccentago deviation, AAf(\og Myir)/Af{log A4^ir) = [Af h s {^og M^ir) — 
AfcPLi^og Myir)]/J^Hs{^'^S ^^vir), ^ function of the halo mass for z from 0.2 to 1.4 in steps of 0.4 (green, blue, purple and red lines). 
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Figure 3. Pecentage deviation, AA/'(log Af„ir)/A/'(log M^i^) as a function of the halo mass for z from 0.2 to 1.4 in steps of 0.4 (from top 
left to bottom right) and e from -3.5 to -7.5 in steps of -1.0 (yellow, green, blue, purple, red lines, respectively). 



The mass function (MF) for this parameter set (which 
we wiU refer to in the foUowing as the fiducial model) is 
plotted in the left panel of Fig. [2] for three different redshift 
values. Although not clearly visible, a kink is present in the 
MF at the scale log Mvir ~ 13 as a consequence of the inter- 
polation scheme adopted to reproduce the MF from N - body 
simulations. Actually, such a feature is not an artifact, but 
rather the evidence of the chameleon effect. Indeed, the scale 
at which the kink is present is the same as the one mark- 



ing the onset of the chameleon mechanism which increases 
the abundance of intermediate mass haloes. This can also be 
seen from the right panel where we compare the HS MF with 
the CPL one. For Mvir » Mth, the chamaleon effect fully 
masks the enhancement of the gravitational force introduced 
by the modified gravity potential thus leading to a MF which 
is the same as the CPL one. In the opposite regime, the ef- 
fective field is no more negligible and the gravitational force 
is significantly boosted leading to a marked increase of the 
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abundance of intermediate mass haloes. Since the strength 
of the effective field is a function of z, the threshold chama- 
leon mass will change with the redshift so that the ratio 
between the HS and CPL MPs does depend not only on the 
mass, but on the redshit too. 

Although discriminating between HS and CPL is inter- 
esting in its own, it is also instructive to look at how the HS 
MP depend on the ,f{R) parameters. This can be inferred 
from Fig.O where we plot the percentage deviation from the 
fiducial model varying e for the fiducial n at different red- 
shifts. As an encouraging result, we find that AAf/M can 
take quite large values in the range 12 ^ logM„ir ^ 15 
for all z. Similarly large values can be achieved for larger 
masses, but we warn the reader not to overcome the re- 
sults in this range since they refer to MPs which have al- 
most vanished so that even a small difference in value gives 
rise to a large AAf /Af. Explaining the shape of AM/M is 
not a straightforward task. On one hand, the scaling of the 
chameleon threshold mass Mth as l/flo - = dex(3e/2) 
implies that deviations from the fiducial MP are positive 
or negative depending on the mass being larger or smaller 
than the corresponding chamaleon mass. On the other hand, 
since Geff also d epends on e, the mass v ariance a{Mvir), 
which enters the ISheth fc Tormen I {l99^ ) profile through 
V = 5c/cr, will be a function of e thus introducing a second 
way through which e affects the MP. ft is actually difficult 
to disentangle the two effects so that the final non mono- 
tonic behaviour in Pig. [3] is not immediate to quantitatively 
explain. It is, on the contrary, much easier to understand 
why we find that n has only a very minor role in setting 
the shape and the amplitude of the MP with AM /N being 
smaller than 1%. Indeed, since the background expansion 
is almost the same f or each n (see, e.g., the discussion in 
ICardone et aTl 120121 ') . the only way n can impact the MP 
is through the effective gravitational constant Geff - Por the 
HS models we are considering, this dependence is actually 
quite mild and is further smoothed out by the procedure 
leading to the final MP. Note that this is partially a con- 
sequence of how we set the fiducial parameters. Indeed, for 
(n, e) = (1.5, —6.0), we get Ci/ci ~ 0.004 so that, in the high 
curvature regime, we get f{R) ~ R — m^ci/c2 with n playing 
a marginal role in setting the ratio Ci/c2. Should we have 
chosen, e.g., (n, e) = (1.5,-4.0) as fiducial, we would have 
obtained ci/cl — 0.4 and hence the term m^(ci/c|)(m^/_R)" 
in the high curvature Lagrangian would have played a much 
significant role likely leading to a stronger dependence of 
the MP on n. However, such a large e does not lead to a 
background expansion in agreement with our fiducial CPL 
one so that we have not considered this kind of models. 

As a final remark, we stress that both Pigs. [2] and [3] 
refer to the theoretical MP. This is not what is actually 
observed so that the differences between the HS and CPL 
MPs or among the different HS model parameter sets can be 
strongly suppressed when the selection effects introduced by 
the particular method adopted to observationally determine 
the MP are taken into account. Should the observational MP 
be defined only for \og M^r > 14 or be measured with too 
large uncertainties at higher z, the possibility to discriminate 
among the HS and CPL models and /or to constrain the HS 
parameters (n, e) would be seriously compromised. 



3 WEAK LENSING PEAKS 

Being the largest and most massive mass concentrations, 
galaxy clusters are ideal sites to look for gravitational lens- 
ing effects. Should the source be aligned with the cluster, 
the formation of spectacular arcs takes place, but, for most 
cases, the main effect is a coherent distortion of the shape of 
background galaxies, which appears stretched tangentially 
around the cluster. The distortion field can be used to re- 
construct the shear map and then infer constraints on the 
cluster mass distribution. On the other hand, a shear map of 
the sky presents distinct peaks corresponding to the cluster 
positions thus offering a technique based on weak lensing to 
identify galaxy clusters. 

As a p eak finder, we c onsider here the aperture mass 
defined by (|Schneider II 19961 ) 



Map(6») = J ti{e)u{i}-e)(fe = j ^t{e)Q{^ -e)d^e (19) 

where k{9) and 7t(6) = — 7?.[7(^) exp (— 2i(j!>)] are the 
convergence and the tangential shear at position 6 = 
(i9 cos (j), ■& sin (j)) and U [■&), Qij}) are compensated filter func- 
tions related to each other by the integral equation 



g(i9) = -C/(i9) + ^ j'^ U{d')'&'dd' . 

In order to detect a cluster as a peak in the aperture mass 
map, we have to preliminarily estimate the Map variance and 
then set a cut on the signal -to -noise (S/N) ratio. To this 
end, we have to specify how we compute both the signal and 
the noise which we address in the two following subsections. 



3.1 Halo model 

Por given lens and source redshift [zi, Zs), the aperture mass 
depends on the mass density profile of the halo acting as 
a lens. We will model c luster haloes using t he spherically 
symmetric NPW profile (jNavarro et al. |[l997l ') : 



p{r) 



{M^ir/^1lRtir)fNFw{c^ir) 



with X = r/Rs, y = r/R^ir and 

/jVFw( 



In ^1 -|- Cvir'} Cuir/(l H~ C-uir) 



(20) 



(21) 



Note that we have reparameterized the model in terms of the 
virial mass M^ir, defined as the mass within the virial ra- 
dius Rvir where the mean density equals AvirPcrit{zi) with 
pcrit the critical density at the lens redshift, and the concen- 
tration Cvir = Rvir/Rs- Accordiug to N-body simulations, 
the NPW model can be reduced to a one parameter class 
since c^ir correlates with the virial mass M-uir- Actually, the 
slope, the scatter and the redshift evolution of the c„ir - M^ir 
relation are still matter of controversy with different results 
available in the literatu r e. As a fiducial case, here we follow 
iMufioz - Cuartas et al. I (|201 j ) setting 
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log C„i: 



(0.0292; - 0.097) log (h-'^A'hir) 



110.001 



+ 



2469.720 



zi + 16.8885 {zi + 16 



(22) 



referring to this in the following as the MCll relation. How- 
ever, in order to investigate the dependence of the results 
on the adopted mass - concentration law, we will also con- 
sider the Cvir - Mvir relation empirically found by Buote et 
al. (2007, hereafter B07) from observed X-ray galaxy sys- 
tems spanning the mass range (0.06, 20) X 10" Mq, namely : 



Cvi- 



Co 



1 + zi VlOi^Mg 



(23) 



with {co,a) — (9.0,-0.172) as inferred from the fit to the 
full sample. It is worth noting that both Cvir - Mvir relations 
are affected by a significant scatter. In order to take this 
into account, one should convolve the aperture mass for a 
given Mvir value with a lognormal distribution centred on 
the predicted Cvir value and with variance equal to the scat- 
ter itself. However, we follow here the common practice of 
neglecting this issue and simply insert the anal ytical expres- 
sion f or the convergence k of the NFW profile (|Bartelmann I 
1 19961 : 1 Wright fc Brainer"d1l2000^ into Eq.lUS} with Cvir set 
according to Eq.((22| or (p3)) . 

A caveat is in order here. Both the NFW profile and the 
Cvir - M-air relations have been inferred from N - body simu- 
lations carried on in a ACDM cosmological framework hence 
implicitly assuming that, on galaxy scales, the gravitational 
potential is Newtonian. On the contrary, we are here investi- 
gating f[R) theories so that, strictly speaking, the classical 
Newtonian theory does not hold anymore and the poten- 
tial is modified. Actually, a Yukawa -like term is added to 
the Newtonian potential with a scale length which, due to 
the chameleon effect, depends on the environment and the 
strength of t he effective field related to the modified gravity 
Lagrangian (| Faulkner et al. I l2007l : ICardone fc Capozziello I 
I2OIII ). However, for the typical values of the f{R) model 
parameters we are going to consider, the effective potential 
is almost identical to the Newtonian one so that we do not 
expect deviations of the halo profile from the NFW one. 
This is indeed what I Schmidt et al. I (|2009l ) have found in 
their analysis of the halo profiles from N - body simulations 
carried on for the same HS model we are considering here. 
Although a Cvir - Mvir relation was not reported, we expect 
that the (unknown) actual relation is within the range set 
by the MCll and B07 relations. 



3.2 Filter function and S/N ratio 

It is worth noting that, being sensitive to all the matter 
along the line of sight, the observed Map is actually the sum 
of the contribution due to both the cluster and the uncorre- 
lated large scale structure projected along the same line of 
sight, i.e. Map = M^i"" -|- M£f ^. In the usual approach, one 

LSS 



^ap ^^-^ap 

considers that, being a density contrast, Map'^ averages to 
zero so that no bias is introduced in the Map statistics. On 
the contrary, the LSS will contribute to the variance repre- 
senting an additional source of noise (iHoekstralbOoH ). The 
filter functional form and its parameters are then chosen 



based on a comparison with sim ulated data and the charac - 
teristics of the survey (see, e.g. iHetterscheidt et al. Il2005l ). 
Actually, such a method is not fully efficient in filtering out 
the LSS contribution and is, moreover, related to the under- 
lying cosmology adopted in the reference simulation. 

A possible way out of this problem relies on the use 
of the so called optimal filter introduced by Maturi et al. 
(2005, hereafter M05) to take explicitly into account both 
the shape of the halo signal and the underlying cosmology. 
Following M05, we set the Fourier transform of the filter as 



(27r) 



m 

Pn{£) 



(24) 



where f is the Fourier transform of the signal (in our case, 
the tangential shear component) and Pn{£) the noise power 
spectrum as a function of the angular wavenumber £. This 
latter is made up by the sum of two terms : 



with 



Pn{£) = Pe+P-,{£) , 



2 n„ 



(25) 



(26) 



the noise contribution from the finite number of galaxies 
(with number density Ug) and their intrinsic ellipticities {a^ 
being the variance), and Pj(£) = ^Pk{£) the noise due to the 
LSS. Note that this latter is set equal to half the convergence 
power spectrum since we use only one component of the 
shear. Under the Limber flat sky approximation, it 



p^w 



V 2c2 



with 



X(2) 



' rxh 
Jo 



c 



Ps -, 
X 



dz' 



dx (27) 



(28) 



the comoving distance to redshift z (with Xh the distance 
to the last scattering surface), VV(x) is the lensing weight 
fucntion 



W(x) = 



1 7 ) Pxix')x'dx' 



(29) 



and Px{x)d-X = Pz(z)d z is the source reds hift distribution 
here parameterized by ([Marian et al. II2OIII ) 



^ Ea.l l27l l is the same as in GR based dark energy models which 
could be surprising at first sight since we are using a modified 
gravity theory. Actually, it has been shown llTsuiikawa I [20o3) 
that, for f{R) theories, the convergence power spectrum is still 
given by Eq.l l27ll provided the matter power spectrum Ps{£/Xy ^) 
is computed taking care of the running eff'ective gravitational con- 
stant and the result is multiplied by the factor l/f'{R). Since, for 
the models we are considering, /' (R) ~ 1 with great accuracy, we 
have neglected this correction when computing P^. 
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Pz[z) = — — 
20 \Z0 



exp 



(30) 



and normalized to unity. Following iMarian et aTl (|201lh . 
we set (/3,zo) = (1.5,0.6) so that Zm = 0.9 is the median 
redshift of the sources as expected for the Euclid mission. 

F ollowing lB"artelman"rn (|l996l ) and I Wright fc Brainerd I 
l|2000l) for the shear profile of the NFW model, one finally 
gets for the Fourier transform of the filtei0 : 



1 


[(- 


(27r)3 









(31) 



where 'Ticrit ~ Ds / {4,-RGDdDda) is the critical density 
for lensing (depending on the lens and source redshift), 
g{cvir) = fNFw{cvir)/cvir, 6a is the augular scale corre- 
sponding to the characteristics radius Rs, and we have de- 
fined the dimensionless quantities 



Jo 



Pn{£) 



■Ml 



(32) 



(33) 



with 7(^) the dimensio nless shear profile given by 
l| Wright fc Brainerd llioool '). Taking the back Fourier trans- 
form of Ea. (|31|l and setting Q = "i! in the aperture mass 
definition, it is only a matter of algebra to finally get : 



1 Map{l^,es) 



with (5 = ■ff/ds and = e/Os) 



(34) 



the halo properties. This is, of course, not the case. First, the 
halo virial mass M^ir and redshift zi directly enter the def- 
inition of 9s which also depends on the adopted c^ir - M^ir 
relation. Second, the halo model sets the dimensionless shear 
profile 7(5) which enters both in the filter derivation and 
the Map function. On the contrary, the dependence on the 
source redshift Zs has been integrated out when computing 
the filter function. This is a consequence of how the optimal 
filter has been constructed, i.e., imposing that the filter is 
optimized to find halos with a given mass profile, virial mass 
and redsift notwithstanding the source redshift. 

This latter quantity, h owever, plays a role when com- 
puting the noise given bjQ (|Maturi et al. |[2005l ) 



1 

2^ 



Pe|*(£)| ed£ 

^crit 

/o [n + (l/2)P.(^)]2 
so that the S/N ratio readfl 



£M (36) 



5(^. p) = _L Mv,r/iTTRl,, 2-Kelg{c^ir) Map{^) ^^^^ 



/2tv 



where p summarizes the parameters the S/N ratio depend 
on, namely the halo quantities {M^ir, zi) and the source red- 
shift Za (entering because of T^aru)- Since we are interested 
in the total S/N ratio, we integrate over the source redshift 
distribution thus finally obtaining 



S{'&;Zl,h'hir) 



S{-&; Zl,Zs, lvhrr)Pz{Zs)dZs 



(38) 



which, for a given Cvir-M^ir relation and mass profile, de- 
pends only on the halo virial mass and redshift. 



Map = / liO^'d^' (35) 
Jo 

^ * [Os + i''^ - 2CC' cos e) ^''^j cos (2e)de . 

A conceptual remark is in order here. Eas. (|34p has 
been obtained starting from Ea. (|19p and assuming 
that the shear profile is the same as the cluster one. 
Actually, the tangential shear in Eq.(|19p is obtained 
from the observed galaxies ellipticities so that it is 
not strictly equal to the model shear which one es- 
timates from the cluster model. However, since we 
are interested in making analytic predictions, such 
an identification is unavoidable. It is also worth not- 
ing that what is actually observable is the reduced 
shear 7/(1 — /t) which is, however, equal to 7 in the 
weak lensing limit (k << 1) we are interested in here. 

A naive look at Eas. p4|l - (35) could lead to the surpris- 
ing conclusion that the aperture mass does not depend on 

^ Note that the minus sign comes out from our conven- 
tion on the sign of the tangential shear component. 



3.3 Cluster detectability 

Eq. (|38|l enables us to estimate the S/N ratio for a NFW 
halo with virial mass M^tr and concentration c„ir acting as 
a lens at redshift zi. In order to compute S, we need first 
to assign the survey characteristics (determining the noise 



^ Note that it is customary to use Pn{£) instead of in 
Ea.([36)l for the total variance to take into account the 
contribution of both the shot noise and the large scale 
structure. However, we are here considering the signal 
as the sum of the cluster and LSS peaks so that only 
has to be considered as noise term. 

A cautionary remark is in order here. The {M^ir,c^ir) 
parameters entering the S/N ratio refer to the cluster 
halo responsible of the signal one is detecting. This halo 
can also be different from the template model used to 
generate the filter. However, the (M^i^^c^ir) template 
parameters setting the filter normalization cancels out 
when computing the S/N ratio so that, in order to sim- 
plify the notation, we have not used different symbols to 
differentiate the template and actual halo parameters. 
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Figure 4. Left. Signal - to - noise ratio as function of the virial mass for different redshift values, namely z = 0.5, . . . , 1.3 in steps of 0.2 
(from the upper blue to the lower black curves) adopting the fiducial cosmological model and the MCll c^ir-M^ir relation. Centre, 
Same as left panel but for the B07 c^i^ - M^ir relation. Right. S/N percentage difference, A{S/N) / (S/N) = (<Smcii ~ Sbo7)/Smcii^ 
as a function of mass and redshift (same color scheme as before) for the fiducial cosmological model. 



power spectrum P^) and set the filter scale i} . We con- 
sider an Euclid - like survej0 (|Laureiis et al. |[201ll ) and take 
a survey area of 15000 deg^ with an ellipticity dispersion 
CTe = 0.3 and a number of galaxies Ug = 30 gal/arcmin^. 
Note that these are the survey goal, but one can easily scale 
the S/N ratio noting that S cx y/n^, while the total number 
of peaks scales linearly with the survey area. 

The choice of the filter scale i? asks for some caution. 
The optimal filter has been designed taking into account 
the NFW profile in order to maximize the signal. A natural 
scale would therefore be i9 = 6*^ since most of the mass 
contributing to the lensing signal is contained within this 
aperture. For a cluster with Mvir = 10^^ Mq at the median 
survey redshift zi = 0.9, we get 6s — I arcmin, while the 
virial radius subtends an angle O^ir — 3.5arcmin (in our 
adopted fiducial cosmology). However, not all the clusters 
have the same mass and the same redshift. Setting = 
1 arcmin would be an optimal choice for these median values, 
but will strongly underestimate the signal for clusters at 
lower redshifts. On the contrary, a varying ■£) would allow to 
maximize the S/N at every redshift, but would introduce an 
inhomogeneity in the computation of the number of peaks. 
As a compromise, we therefore set i9 = 2 arcmin noting 
that the filter actually cuts the contribution from regions 
outside few times 63 and vanishes on scales larger than 9vir 
(which can be smaller than 2 arcmin for low mass and/or 
high redshift clusters). 

It is instructive to look at the S/N ratio as a function of 
the cluster redshift and virial mass for a given cosmological 
model. To this end, we consider the fiducial HS model with 
(nj\/, n, e) — (0.273, 1.5, —6.0) and look at how <S depends on 
the adopted c^ir - Mvir relation. Fig.|4]shows S as function of 
the virial mass at different redshift values for the MCll and 
B07 cases (left and central panels, respectively). For both 
models, the trend with {z, log Mvir) is the same : the higher 
the redshift, the larger must the virial mass be in order to 
attain a given S value. This can be qualitatively explained 
by noting that the integration in Ea. (|38|) is performed on a 
smaller interval with increasing z thus reducing 5. In other 
words, the larger is the cluster redshift, the lower is number 
of sources available for lensing so that the S/N ratio takes 
smaller values. In order to compensate for this reduction. 



^ |http : //wmi . euclid-ec ■ org | 



one has to increase the virial mass since S approximately 
scales as M^ir with the slope a depending on the adopted 
Cvir - Mvir relation. 

The right panel in Fig. |4] compares 5 for the MCll and 
BOl cases as a function of Mvir and Zi showing that the 
S/N ratio turns out to depend strongly on the Cvir-Mvir 
relation. In particular, for less massive systems, the MCll 
relation predicts a values of 5 which can be up to 100% 
larger than the B07 ones. The difference is, however, not 
monotonic reverting its sign for intermediate mass haloes 
and becoming quite small for very large mass. Such a strong 
sensibility to the Cvir - Mvir model is not unexpected and is 
a consequence of the different concentrations for the same 
virial mass. In particular, the ratio CviJ /c^ir^^ depends on 
both zi and Mvir thus explaining why the A5/5 may change 
sign according to which model predicts the larger Cvir value. 

As a final remark, we stress that the above results are 
almost independent of the underlying cosmological model. 
This can be easily explained by looking at how cosmology 
enters in the evaluation of the S/N ratio. First, we note that 
S depends on the integrated Hubble rate through the lensing 
critical density T^crit- This quantity only weakly depends on 
cosmology since it involves an integral of a ratio of distances 
over the source redshift. On the other hand, the convergence 
power spectrum Pk{£) enters the filter definition and hence 
both A4ap and aap so that its impact is reduced when consid- 
ering the ratio of these two quantities. Moreover, for the HS 
model parameter space we are interested in, the Hubble rate 
H{z) is almost the same as the CPL one so that changing 
the {Q,M,n,£) values has a very minor impact on the S/N 
ratio. This is indeed what we have checked. Varying (f2j\/, e) 
shifts S by less than 10% (but typically a shift larger than 
1% is difficult to get), while the S/N ratio may be considered 
independent of n within a very good approximation. 



4 PEAK NUMBER COUNTS 

Having detailed how the S/N ratio can be computed and 
having determined the theoretical MF, we can now estimate 
the number density of haloes with mass Mvir that produces 
significant peaks in the aperture mass map. To this end, we 
first have to take into account that a S/N threshold for the 
weak lensing signal does not correspond to an equally sharp 
threshold in halo mass because of the scatter in A4ap caused 
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Figure 5. Left. Total number of detectable peaks as a function of the threshold S/N ratio for the fiducial HS model adopting the MCll 
Ciiir - ^Ivir relation; solid, dashed and dotted lines refer to the total number and the contribution from clusters and LSS, respectively. 
Centre. Same as before but showing the total number only (in the regime dominated by the clusters contribution) for the MCll (red) 
and B07 (blue) mass - concentration relations. Right. Percentage deviation, ANpi^/Npi^ = i^^pff — ^pk'^)/^pif ''^ function of Sth 
with the labels HS and CPL referring to the fiducial HS and CPL models. 



by the shot noise from discrete background galaxy positions 
and the intrinsic ellipticity distribution. A halo of mass M^ir 
has therefore a certain probability p{Map\Mvir) to produce 
an aperture mass Map which we can model as a Gaussian : 



p{Map\Mvir) oc exp 



Map - MapiMyir) 



(39) 



where Map{Mvir) is the theoretically expected value and 
ffap the variance in Ea. (|36|l integrated over the source 
distribution. The probability that the 5'/A'^ ratio will be 
larger than a given thre shold value will then be given by 
l|Bartelmann et al. llioO^ ') : 



p{S > Sth\Myir,z) = 2*^''^'^ 



SiM^ir 



^/2 



(40) 



so that the number density of haloes giving a detectable 
weak lensing peak finally reads 

MensiMvir, z) = p{S > Sth\Mvir, z)M {M^ir , z) (41) 

where, hereafter, we will drop the label I from z to denote 
the cluster redshift. The total number of peaks generated by 
cluster haloes and with S/N larger than a threshold value 
Sth is then obtained by integrating over z and multiplying 
for the survey area Q. thus finally reading : 

r^(z) f°° 
E(z) Jo 

where we set {zl,zu) = (0.1,1.4) as redshift limit^, while 
r{z) — {c/Ho)^^x{^) is the comoving distance. 

The number of observed peaks is, however, larger than 
J^haio because of the contamin ation from the LSS. This lat- 
ter term may be estimated as (|Maturi et al. II2OI0I . [2OI1I 1 



° Note that the survey will likely detect galaxies over a much 
larger range, but we have cut the redshift range since, as will be 
shown later, the number of peaks outside this range is negligible. 



Ml 



' (27r)3/2 V -ap 

with Kth = SthCTap and 



■ exp 



1 f Hth 

2 V CTau 



(43) 







2 , 

fdl 




4>{£) 


^ £d£ 



(44) 



Note that A/lss only depends on the noise properties and 
the threshold S/N ratio, but not on the lens mass and red- 
shift since it is not related to the particular cluster at hand, 
but to the LSS. For this same reason, A/lss is determined 
by the matter power spectrum (and hence the underlying 
cosmological scenario) entering Pn{£). This will actually of- 
fers a further opportunity to discriminate among GR based 
dark energy models and modified gravity theories. The total 
number of observed peaks will finally be given by : 

jVpk{S > Sth) = JVhaloiS > Sth) +JVLSs{Sth) ■ 

A look at Fig.[5]helps us to highlight many important issues. 
First, in the left panel, we plot Npk{S > Sth) as a function 
of the threshold S/N ratio for the fiducial HS model and 
the MCll Cvir - Muir relation showing separately the contri- 
bution of the clusters and LSS terms. As expected, Npk is 
dominated by the LSS term for small 5'/A'' ratio, that is to 
say, the smaller is the S/N ratio, the larger is the probability 
that the detected peak is a fake due to the LSS rather than 
the evidence for a cluster. This is in agreement with com- 
mon sense expectation and previous analysis in literature 
using different cosmologica l models and survey parameters 
(|Hetterscheidt et al. I [20051 : iMaturi et al~ll2010l '). The term 
A/lss, however, quickly becomes subdominant so that one 
can confidently be sure that all peaks detected with a S/N 
ratio larger than Sth — 5 are due to clusters. Note that we 
find a va lue for Sth comparab le but larger than what is sug- 
gested in lMaturi et al. I (|201G| 1 because of differences in both 
the cosmological model and the survey characteristics. 

The central panel in Fig. [5] shows the impact of the 
adopted c^ir - M^tr relation zooming on the high S region 
which is dominated by the clusters term. It turns out that 
the number of peaks is systematically larger if one adopts 
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Figure 6. Left, Total number of peaks in redshift bins for the HS (solid) and CPL (dashed) models adopting the MCll mass- 
concentration relation and threshold S/N ratio, Sth = 5,7,9 (red, purple, black lines). Right. Percentage difference as a function of 
2 for the three Sth values adopted. 
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Figure 7. Percentage deviation of Npi^{z) as a function of z with respect to the fiducial HS model with the MCll c^ir - M^ij. relation. 
In each panel, only one of the parameters (n, e, fjjvf , erg) is changed. Top left. Results for e = —4.0 (green), e = —5.0 (blue), e = —7.0 
(purple). Top right. Results for n = 1.0 (blue) and n = 2.0 (purple). Bottom left. Results for Hjv/ = 0.234 (green), = 0.254 (blue), 
Hm = 0.294 (purple). Bottom right. Results for trg = 0.61 (green), erg = 0.71 (blue), erg = 0.91 (purple). 



the B07 rather than the MCll relation. This is actually not 
a fully unexpcted result. Indeed, as the right panel of Fig.|3] 
shows, the S/N ratio is larger for the B07 model in the high 
mass regime. Put in other words, for a given Sth, the num- 
ber of high mass clusters with 5 > Sth is larger for the B07 
Cvir - Mvir relation. Since these are the main responsible of 
the total number of high S/N peaks, it is then straightfor- 



ward to explain why we indeed find that Npk is larger for 
the B07 relation and why the offset increases with S. 

The most interesting issue to address is whether the 
number count of detectable peaks in the aperture mass maps 
allows to discriminate between the HS and CPL models 
notwithstanding their equal expansion rate. The right panel 
of Fig.[5]offers a first answer to this question. S hould we rely 
on clusters peaks setting 5 > 5 for detection (jCeller et al. I 
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I2OIOI : iKurtz et al.1 \20li ). the difTerence between the two 
models amounts to a modest 3% with the HS case giv- 
ing a larger number of peaks. This is somewhat surpris- 
ingly considering that the enhanced growth of structures 
in f{R) theories should lead to a larger number of clusters 
for a given mass. Actually, such a result is not fully unex- 
pected if one looks back at the right panel in Fig. (2] Here, we 
have shown that the HS MF is larger than the CPL one for 
M-uir < 10^* Mq, while the difference becomes negligible at 
higher masses because of the onset of the chameleon effect. 
Since the S/N ratio for haloes with masses below 10^"* M0 
is quite small (unless they are at low z), their contribute to 
the total number of peaks is strongly suppressed thus lower- 
ing the difference with the offset between the HS and CPL 
predictions. It is, however, worth noticing that this same off- 
set significantly increases (up to ~ 15%) and changes sign if 
we decrease the threshold S/N ratio. This is related to the 
onset of the LSS dominated regime. In such a case, what is 
important is the difference in the convergence power spec- 
trum entering the ratio (Jlss /(^ap, while -Pk(^) plays a minor 
role in Miens ■ We can therefore conclude that the number of 
low S/N peaks is a valuable tool to discriminate between 
HS and CPL models. It is also worth stressing that the off- 
set between the two models predictions is (within a good 
approximation) independent of the adopted c„ir - M^ir rela- 
tion as can be seen from the closeness of the red and blu 
curves in the figure. This is an expected result consider- 
ing that this ingredient only enters the determination of the 
S/N , but does not change how the number of peaks depend 
on the underlying cosmological model. 

The total number of peaks is evaluated by integrat- 
ing over the full redshift range thus smoothing out the de- 
pendence of the MF on z. Actually, how the MF changes 
with 2 depends on the underlying cosmological model so 
that it is worth investigating what can be learned by con- 
sidering the number of peaks in redshift bins. This is 
still given by Ea. (|43|) provided one replaces {zl,zu) with 
(2 - A2/2, z A2/2). We set A2 = 0.1 and show the re- 
sulting Npk{z,S > Sth) in Fig. ((6)1. Note that here we have 
turned off the LSS term. Since these peaks are fake detec- 
tions due to the LSS, it is obvious that no redshift can be 
attributed to them. We can therefore assume that they have 
been manually deleted from the sample and simply se10 
Mpk{z) = Nhaio{z)- The left panel suggests that binning 
the data has not improved the constraints, but this is only 
part of the story. Indeed, the percentage difference between 
HS and CPL predictions are quite modest (~ few %) up to 
2 ~ 1, but they become as large as 20% in the highest red- 
shift bins. Such a behaviour can be qualitatively explained 
looking back at the right panel of Fig.[2] showing that the 
larger is 2, the larger is the difference between the HS and 
CPL MF. Probing the high -2 end of the Npk{z) curve is 
therefore an indirect method to investigate the MF in the 
regime where it is most sensible to the underlying gravity 
theory. As a further remark, we also note that the percent- 



Such an assumption is actually also mathematically motivated. 
Indeed, we are now considering dj\fpk/dz. Since J\fpk is the sum 
of the two terms due to the haloes and the LSS, respectively, 
and the LSS one is independent on 2, we trivially get dMpu/dz = 
dMhalo/dz + dMLSs/dz = dMhalo/dz = Mhaloiz). 



age difference between HS and CPL predictions at a given 
2 is smaller for higher S/N thresholds. This can be quali- 
tatively explained noting that, in order to have a large S 
peak, the cluster must be very massive so that the larger is 
Sth, the more we are investigating the high mass tail of the 
MF. As already stated before, the HS and CPL MFs become 
degenerate for M-uir > 10^"* Mq so that a sample probing 
this range will be unable to discriminate between the two 
competing models. 

Although discriminating between CPL and HS scenar- 
ios is likely the most intriguing issue to consider, it is never- 
theless of deep interest to investigate how Npk{z) depends on 
the HS model parameters. To this end, we hereafter adopt 
the MCll relation since this gives a lower number of de- 
tections thus providing conservative estimates . Fig.[7]shows 
that only e has an important effect on the peak number 
with deviations which can be of order 10 — 30% depend- 
ing on the redshift bin considered. Comparing Fig. [7] with 
Fig. [3] clearly shows the impact of moving from a theoretical 
to an observable quantity. Indeed, the deviations from the 
fiducial model are almost halved the expected values, be- 
cause only the haloes with masses large enough to lead to 
a S/N ratio larger than the threshold can now be detected 
thus removing the part of Fig.|3]with large deviations. As a 
consequence, while the possibility of discriminating among 
different e values is still significant, the dependence on n al- 
most disappears with deviations being smaller than 1% and 
only present in the highest redshift bins. Finally, although 
not f{R) parameters, we also plot the deviations from the 
fiducial model due to changes of the matter density Qm and 
the variance as- The bottom left panel of Fig. [7] shows that 
deviations as high as ~ 40% can be obtained by changing 
Qm, but the way ANpk/Npk. depends on z is different from 
what happens when varying e thus suggesting that no strong 
degeneracy between these two parameters is present. Simi- 
larly, the lower right panel tells us that erg also plays a major 
role in determining the peak number counts in the highest 
redshift bins thus suggesting that strong constraints could 
be put on this parameter. 



5 FISHER MATRIX FORECASTS 

In order to quantify the power of peak number counts to con- 
strain the HS model parameters, we carry on a Fisher ma- 
trix analysis considering as observed data the total number 
of peaks with 5 > Sth in equally spaced redshift bins cen- 
tred on 2 and with width A2 = 0.1 over the range (0.1, 1.4). 
As a first step, we need to define a likelihood function jC 
to quantify how well a theoretical model matches the ob- 
served number counts. Should the errors be Gaussian, one 
can resort to the usual statistics and define —2ln£. — x^- 
However, when dealing with number counts, one can assume 
Poisson erro r jf] and then rely on the so called C statistics 
jCash II 19791) to define : 

* Note that, if there are uncertainties in assigning an object 
to a given bin, one should better resort t o the pull statistics 
llFogli et al. II2OO2I: ICampanelli et al. Il201lh . If we assume that 
the peak redshift has been determined with sufficient accuracy, 
we can neglect this possibility and simply rely on the usual ap- 
proach. 



Peaks statistics and f{R) theories 13 



0.91 ■ 



b 0.81 ■ 



0.71 ■ 



98.5 
?J 1.5 
-97.5 



-3.32 
^ -6 
-8.68 



— 1 1 r- 



0.263 




H 1 1- 



H 1 H 





H 1 1- 




— 1 1 r- 




—I I i_ 



-97.5 1.5 98.5 
n 



-8.68 -6 -3.32 
e 



— S/N=5 

— S/N=4 



0.273 



0.283 



Figure 8. Fisher matrix constrain in 2D spaces marginalizing over all the parameters but the ones on the axes. 
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where p is the set of model parameters, A/tin the number 
of redshift bins and, for notational clarity, we have defined 
Xi = 7Vpfc(zi,p) and m = N°ji''{zi) to denote the theoret- 
ical and observed number of peaks in the bin centred on 
Zi (and with width Az — 0.1). The Fisher matrix elements 
will then be given by the second derivatives of the logarithm 
of the likelihood with respect to the parameters of interest 
evaluated at the fiduc ial values. Starting from (|45[l . one gets 
jHolder et al. Il2001al ') : 



E 



dp. 



(46) 



where A^''* is the expected number of peaks in the A: - th bin 
for the fiducial model. The inverse of the Fisher matrix gives 



us the covariance matrix with its diagonal elements repre- 
senting the lowest variance that one can achieve on the pa- 
rameter Pi. It is worth noting that the Fisher matrix approx- 
imates confidence regions as a Gaussian ellipsoids, while the 
true confidence o nes can have b r oad tai ls or significant cur- 
vature. However, [Holder et all (|2001al ) have checked that 
this is not the case for number counts by comparing with 
Monte Carlo analysis of simulated datasets. We expect that 
our estimated iso - likelihood contours provide a reliable ap- 
proximation of the degeneracies in the parameters space. 

Since here we are mainly interested in how peak num- 
ber counts constrain f{R) theories, we will not consider the 
full eight dimensional parameter space. On the contrary, we 
will set p = {^Im , e,n, as), while {h,WQ,Wa,nps) are held 
fixed to their fiducial values. While this is not what one will 
actually do when dealing with real data, it is worth noting 
that here we are only considering peaks statistics as a con- 
straint. Actually, the background quantities {h,wo,Wa) will 
be severely constrained by future SNela data, while nps is 
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Table 1. Fisher matrix forecasts for the accuracy on the pa- 
rameters (f2i\/, e, n, (Tg) from fitting the Npk{z,S > Sth data for 
different threshold S/N ratio Sth- 





Qm 


e 


n 


o"8 


2 


0.0051 


0.0452 


1.12 


0.068 


3 


0.0007 


0.2523 


6.33 


0.033 


4 


0.0009 


0.8189 


18.8 


0.042 


5 


0.0020 


1.0738 


38.8 


0.034 



determined by the fit to the CMBR spectrum. In a fully re- 
alistic approach, we should therefore both explore the full 
8D parameter space and add the likelihood terms related to 
SNela and CMBR data. In order to better highlight the im- 
pact of peak statistics, we prefer to deal only with this kind 
of data and limit the parameter space. We therefore only 
consider (Hm , e,n, as) since they are intimately related to 
the mass function and the f{R) model. The errors on these 
four quantities turn out to depend on the threshold S/N 
used and are summarized in Table[T] while Fig.[5] shows the 
iso - likelihood contours in 2D spaces marginalizing over all 
the remaining parameters for the cases Sth = (4, 5). 

It is interesting to note how the constraints scale with 
the threshold S/N. On the one hand, increasing Sth has the 
effect of reducing the overall number of peaks thus one ex- 
pects a worsening of the constraints because of the poorer 
statistics. This is indeed the case for (e,n), but not for 
(Hm, eg). This is likely related to a second effect of a larger 
Sth- As a consequence of the mass and redshift scaling of 
the S/N ratio shown in Fig.|4l varying Sth implicitly leads 
to probing a different MF regime. Since the dependence of 
the MF N{Mvir,z) on (yiM^ag,) is different depending on 
the (M„ir, z) region of the parameter space one is consider- 
ing, the constraints will turn out to depend in a complicated 
way from which kind of peaks one is detecting. 

Our main interest in the peaks statistics was moti- 
vated by constraining the ,f{R) model parameters. Taken 
at face values, the constraints on e could look not so en- 
couraging. We stress, however, that they are still competi- 
tive with those from an analysis based on background evo- 
lutionary probe s as c an be appreciated remembering that 

1 Cardone et al. I (|2012h obtained a{s)/e ~ 35% from fitting 
SNela -I- GRBs -I- BAO -I- CMB distance priors vs the present 
forecasted a{e)/e ~ 16% value (for Sth ~ 5) using weak 
lensing peaks counts only. 

As it is apparent from both the numbers and the plots, 
the only parameter which remains fully unconstrained is the 
slope n of the HS model. We could naively anticipate this re- 
sult by simply looking at the top right panel of Fig.[7] where 
one can see that Npk{z) changes by less than 0.1% up to 

2 ~ 1 so that this parameter has a negligible impact on the 
observed peak number. This conclusion can be qualitatively 
explained by looking at the effective gravitational constant 
Qeff in EQ. p5|) . For the e values typical of viable HS mod- 
els, the dependence of Qeff on n almost cancels out from the 
ratio so that the power spectrum Ps{k, z) is almost indepen- 
dent of n with the residual sensitivity being smoothed out 
from the integration needed to get the variance a^{Mvir). 
It is therefore not surprising that the Fisher matrix analysis 



finally leads to a large uncertainty on n. We, nevertheless, 
caution the reader that this is partly a consequence of how 
we have set the fiducial model. Indeed, as already hinted at 
in Sect. 2.3, for (n, e) = (1.5,-6.0), the f{R) Lagrangian is 
quite similar to the ACDM one with the term depending on 
n giving an almost negligible contribution. Investigating how 
the Fisher matrix forecasts change with the fiducial model 
is, however, outside our aims here. 

The stronger dependence of Npk{z) on e which can 
be read from the top left panel of Fig. [3 on the contrary, 
helps us to understand why peaks statistics is so efficient 
in constraining this parameter. Again, we can look at Geff 
to qualitatively explain why the results so strongly depend 
on e. First, we note that, for the values of interest, we get 
Qeff — G/{1 — 10^) within a good approximation. Second, 
the interpolation between gfoG and (jgr in Ea. (|17|) is medi- 
ated by a mass scale which depends on e. Small variation in 
this quantity can therefore significantly change the number 
of peaks thus leading to stronger constraints. 

It is worth stressing that the constraints in Table[T]are 
actually quite optimistic since they are based on the as- 
sumption that one is able to partition the detected peaks in 
redshift bins. First, we note that, should the peak be due to 
the LSS, it will not have any actual counterpart in optical 
images so that the notion of redshift can not be applied. We 
are therefore forced to consider only the results for the case 
Sth = 5 (or larger). Even limiting our attention to these 
high S peaks only, there is still the problem of how to di- 
vide them in redshift bi ns. Ideally , one should rely on optical 
clusters finders (Postman et al. ([l996; K ocstcr ct al. 2 00^ ; 
iMilkeraitis et al. liioiol : iBellagamba et al. II2OI1I ') to match 



the peaks catalog with optical detections and then infer the 
redshift. To this end, it is worth noting that most 
of the optical finders techniques rely on the same 
kind of data as those available for Euclid so that the 
match with weak lensing detections could be imme- 
diate. Moreover, since the cluster redshift estimate 
is based on more galaxies at once supplemented by 
geometrical constraints, the precision on z will be 
similar (if not higher) than that of the photometric 
redshift of a single galaxy. In a first approximation, we 
can therefore assume that the optical finder provides a prob- 
ability distribution function for z which can be modelled as 
a Gaussian with negligible bias and variance ctz = ao{l + z). 
In order a peak with Zi ^ z ^ Zi+i not to be incorrectly as- 
signed to a different redshift bin, we can roughly ask that the 
3cr uncertainty on z is smaller than the bin width. For our 
assumed value Az — 0.1, this translates in ^ 0.03(1 + z). 
Such a precision could be likely achieved if z is spectro- 
scopically measured (cr^ < 0.001(1 -|- z) according to the 
Euclid red book), but could be too demanding if one re- 
lies on photometric redshift methods {a^ < 0.05(1 -I- z) for 
Euclid). In this second case, one should repeat the above 
Fisher matrix analysis adding a non Poissonian uncertainty 
on Npkiz,S > Sth) and resorting to t he pull statistics 
(|Fogli et al. ll2002l : ICampaneUi et al. 1l20llh . Whatever is the 
method adopted for inferring the peak redshift, it is likely 
that its precision will also depend on the peak S/N ratio so 
that the net effect could be included in our analysis by con- 
volving the theoretically computed Npk{z) with an empiri- 
cally determined selection function. Investigating this issue 
is outside our aims here so that we only warn the reader that 
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Figure 9. Percentate deviation of the total number of peaks from 
the fiducial model for e = —4 (green), —5 (blue), —7 (purple). 



the numbers in Tab le[T] should be taken as (Ukely optimistic) 
upper limits. 

As an alternative, one could rely on the total num- 
ber of peaks as only observational constraint. As shown in 
Fig-El this quantity has a significant dependence on e so 
that it can be a valuable help to narrow down the range 
for this parameter. However, a single quantity is unable to 
put constraints in a 4D space so that we have not carried 
out a Fisher matrix analysis for this case. We stress, how- 
ever, that meaningful constraints could be obtained adding 
background probes such as, e.g., S Nela and GRBs . As a n 
example, we can remember that in ICardone etHI (|2012h . 
we found 7.42 < e -3.48 at the 68% CL. Should one de- 
termine Npk{S > Sth) to be consistent with the fiducial Npk 
value within 10% for Sth = 2, one could safely exclude the 
three models in Fig. [5] thus greatly narrowing the confidence 
range. A joint analysis is, however, needed to investigate the 
role of degeneracies with other HS model parameters. 



6 CONCLUSIONS 

The unexpected discovery of the cosmic acceleration 
promptly raised numerous papers proposing different can- 
didates to drive the accelerated expansion. It soon became 
clear that one can also go beyond General Relativity to 
modified gravity theories. In particular, f{R) theories were 
shown to be able to provide exactly the same background 
expansion of dark energy models so that the issue nowa- 
days is, more than checking their validity, discriminating 
between dark energy and modified gravity. Here, we have 
shown that this is in principle possible by relying on the 
peak statistics, i.e. the number of peaks in the weak lensing 
maps constructed using the mass aperture statistics. In or- 
der to quantify this possibility, we have carried on a Fisher 
matrix analysis to estimate the accuracy on the model pa- 
rameters which an Euclid - like survey can achieve using peak 
number counts as the only observational constraint. 

Although the results are quite encouraging, it is worth 
asking whether they can be improved and refined. First, we 
note that the theoretical quantity we have actually relied 
on is the mass function. One can repeat our analysis by us- 
ing different tracers of the mass function itself such as, e.g.. 
X-ray catalogues. Number counts of X-ray selected clus- 
ters have indeed been widely investigated as a possible tool 



for precision cosmology jHolder et al. Il2001bl: iMantz et al~l 



I2OO8I : iBasilakos et al. II2OI0I : IPierre et al. II2OIII ') so that one 
can naively believe that a joint analysis of the mass func- 
tion of X-ray and weak lensing selected clusters can im- 
prove the constraints on the f{R) model parameters. Ac- 
tually, validating such a prediction is not straightforward. 
On the one hand, since both kinds of observations probe 
the same quantity it is indeed possible that degeneracies in 
the model parameter space are not broken so that no signi- 
ficative improvement is achieved. On the other hand, while 
lensing probes the full matter distribution notwithstanding 
its dynamical state. X-ray number counts are subject to 
the uncertainties related to the cluster mass determination 
which typically relies on scaling relations which have been 
determined in a General Relativity framework. Although 
great care has to be used to take into account this prob- 
lem, we nevertheless believe that such a task is worth being 
carried out in order to complement the peak statistics. 

As a second issue, we remind the reader that we have 
only investigated a subset of the full 8D parameter space. 
Allowing for a larger number of parameters to be varied 
impacts the constraints by introducing further degeneracies 
and hence widening the confidence ranges. In order to com- 
pensate for this degradation of the constraining power of the 
method, one can complement peak statistics with different 
dataset. From this point of view, SNela and CMBR are ideal 
tools. Indeed, the SNela Hubble diagram is a reliable tracer 
of the expansion rate over approximately the same redshift 
range covered by the peak number counts; this offers the pos- 
sibility to severely constrain those parameters most related 
to the distance vs redshift relation such as {h, wq, Wa)- Sim- 
ilarly, the CMBR anisotropy spectrum strongly depends on 
(n, (Tg) so that fitting this dataset allows to set these param- 
eters and hence break the {ClM,as) and {crs,e) degeneracies 
shown in Fig.[8l 

As a final comment, we would like to stress that cosmic 
shear tomography can be efficiently added to peak num- 
ber counts since it is has been shown to be particularly ef- 
ficien t at constraining the HS parameters (|Camera et al. I 
I2OIII ). Considering that both the shear power spectrum and 
the peak number counts rely on the same underlying phe- 
nomenon (the lensing distortion of the images of background 
galaxies from an intervening mass distribution) and both 
will be measured from the Euclid mission, we end up with 
the intriguing possibility to use a single mission and a single 
probe to give a definitive answer to the up to now unsolved 
dark energy vs modified gravity controversy. 
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